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Abstract 

Molecular dynamics studies of water and ions around a charged 
BtnaJl compact ^obular object are described. The hypothetical molec¬ 
ular ion is used as the globular model. It is found that the cesium 
cation is in physical contact on the spherical surface much of the time 
with an incomplete hydration shell. The chloride ion spends part of 
the time in a solvent-separated adsorption configuration and part of 
the time in the diffuse r^on. No adsorption of the sodium ion onto 
the spherical surface is observed. 

1 Introduction 

This communication describes the first molecular dynamics study of water 
and ions around a charged small compact globular object about 1 nm in 
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diameter and includes accurate and efficient treatment of all electrostatic 
interactions. The goal of this work is to gain some basic understanding of 
the organization of water and ions in electric double layers around ’spheri¬ 
cal’ objects since this is relevant, albeit on a somewhat larger scale, to the 
microscopic environment around small colloid particles, compact proteins 
and polyelectrolytes. In previous studies the Poisson-Boltzmann equation 
for large spherical ions has been studied analytically by calculating the sin¬ 
glet and pair correlA>tion function using the Omstein-Zemike equation, as for 
example in the « ^2 ^^nderson and coworkers [1, 2, 3, 4]. In the present 

article, we choose eo .acule as the prototype ’structured sphere’ because 
it possesses a comp’'t> orally manageable size and has a characterized set 
of interaction parameters ba^ed on a considerable number of molecular dy¬ 
namics simulations [5, 6, 7, 8, 9]. This paper is not meant to be a study of 
actual Ceo ions in solution. 

So far, computer simulations of electrochemical double layers have been 
mainly performed for aqueous electrol)rte solutiois near flat structureless 
surfaces [10,11,12,13,14,15] and flat atomically struct’ued surfaces [16,17, 
18, 19]. FVom a computational point of view, flat surfaz-cr can be thought to 
be spheres with infinitely large radius so that one anticipates some common 
features between the two types of systems. In particular the first peak in 
the radial distribution function of a large sphere goes over into the density 
probability function in the coordinate z measured along the surface normal. 
Likevnse a contact ion pair goes over into contact adsorption (also called 
spedfic adsorption ) on the flat electrode surface, though in the case of metals 
the electrostatic image interactions differ. Because of the surfsM^e curvature 
changes, one also anticipates some features that differ in each case. 


2 Molecular Dynamics Method 

Constant (JV, V, T) molecular dynzunics simulations are carried out for two 
dilute aqueous solutions. The first is composed of 502 ST2 water molecules 
[20], four sodium cations, two chloride anions, amd one Cso molecule on which 
eadi C atom carries the charge —e/30 so that the total on the molecular ion 
is —2e. The second is composed of 506 ST2 water molecules, two cesium 
cations, and one Ceo molecular ion. For both the cases, the simulation cell 
has the dimensions of 2.5054 X 2.5054 x 2.5054 nm^, and periodic boundary 
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conditions in three orthogonal directions are used. Each simulation lasts 
1 ns, with the first 200 ps being used to equilibrate the system. Detailed 
description of the simulation method was ^ven in previous papers [11, 12, 
13, 14, 15]. 

Before proceeding, we remind the reader that Ceo or buckminsterfullerene 
[21] is a soccer ball-like molecule with 60 carbon atoms sitting on the vertices 
of a truncated icosahedron. E^ach carbon atom interacts with the atoms 
comprising the solution components through Lennard-Jones potentials and a 
CJoulomb potential. The ion-ion, ion-carbon, and water-carbon interactions 
are similarly modeled. The Lennard-Jones parameters of the carbon atom 
are taken from Steele’s work on graphite [22]. The ion-water interactions 
are modeled by the 12-6 Lennard-Jones function between the ion suid the 
oxygen atom of the water and the Coulomb potentials among charged sites. 
The Lennard-Jones parameters of the ions are taken from Heinzinger’s results 
123]. 

3 Adsorption and Hydration 

The adsorption of the ions onto the structured sphere is partially described by 
the radial distribution fimctions (RDFs), guaceo, Qciceo, and gcaceo- 

If there is adsorption onto sfiecific sites like the. hollow site of a six membered 
ring then simple radial distributions based on the assumed spherical symme¬ 
try will not be adequate. However they clearly represent a good starting 
point for initial discussions. As shown in Fig. 1, there is a distinct first hy¬ 
dration shell surrounding the structured sphere with its highest local density 

0.65 nm away from the center of mass of Ceo and the minimum around 
0.74 nm which gives a coordination number of ~ 38 for the system containing 
sodium and chloride ions and ~ 39 for the system containing cesium ions. 
The two RDFs for water are remarkably similar especially when we realize 
that, as will be discussed later, the RDF for Cs+ ion shows one ion to be 
contact adsorbed. 

The ion-Cso radial distribution functions indicate that the Cs*^ ion makes 
a contact p£ur whereas, in contrast, the smaller ions Na’*' and Cl~ do not and 
appear in diffuse regions only. We discuss the cesium cation RDF first. In 
Figure 1 the bottom curve shows that the RDF for Cs"^ consists of a sharp 
peak at approximately 0.7 nm followed by a broad diffuse region devoid of 
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structure. The Cs'*' peak is consistent with a contact adsorbed ion because its 
position relative to the first water peak is equal to half the difference in the 
Lennard-Jones radii for water oxygen (<r = 0.310 nm) and cesium (<r = 0.392 
nm). This strong adsorption tendency is not unexpected because there is a 
Coulombic attraction between the Cs*^ cation and anion and the large 
size of the cation means a weaker hydration energy and unfavorable mixing 
entropy [24]. 

The RDF for chloride ion-Cao is broad and diffuse with the hint of a wide 
mskximumnear 0.85 nm and a weak minimum at approximately 1.05 nm. The 
distance of closest approach is approximately 0.65 nm, indicating that the 
anion occasionally touches the surface of Ceo but does not form a stable con¬ 
tact pair. The broad peak in ^acso lying between the first and second peaks 
in the water-Cm RDF, gtiiOceo, could be expluned by some small tendency to 
take part in solvent-separated adsorption on the sphere. The solvent-induced 
interaction is attractive for the two like ions [25]. If we assume the chloride 
ions, whose distance firom the center of mass of Cso is less th^m 1.05 nm where 
a weak minimum in ^aceo is found, as solvent-separated adsorbing species, 
then the probability for such configurations would be about 26%. 

The radial distribution function for the sodium cation-Ceo in Figure 1 
shows no probability for an approach closer than approximately 0.75 nm. At 
this point the Cso anion and Na'*' cation are separated by one water molecule. 
At the position of the first weak maximum in yNaCso at approximately 0.9 
nm the separation is equivalent to two water molecules. We conclude there¬ 
fore that although the Coulomb interaction between the sodium ion aind the 
Cso ion is attractive, the probability for observing a solvent-separated ad¬ 
sorption of the cation onto the anion is very small, since the cation is fully 
hydrated and the sodium-water interaction is very strong. By comparing 
the adsorption behavior of Na*^ with that of Cl~, it is clear that the size of 
the ion plays a critical role in the adsorption process. Theories and many 
experimental results have shown that larger ionic size corresponds to a more 
negative entropy change and a less negative enthalpy change when the ion is 
transfered from nonpolar solvent to the aqueous environment [26]. In other 
words, large ions behave hydrophobically [27] and have a stronger tendency 
to be excluded from the solution. 

Figure 2 describes the hydration structure of the ions through the ion-0 
and ion-H radial distribution functions. The central pjur of plots for Na-H 
and Na-0 in Figure 2, shows the sodium ion to be most strongly hydrated 
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with distinct first and second shells. The hydration structures of the sodium 
cation and the chloride anion are not significantly perturbed by the Cm 
molecular ion and are very similar to the ions in the bulk NaCl solution 
[28]. On the other hand, the hydration shell of the cesium ion is not ailways 
complete and the average hydration number reduces from 8.3 in the absence 
of the Cm [28] to 7.5 in the presence of the Cm- The values for the first peaks 
of gcM and gcMt 3.4 and 2.2, in Figure 2 are much smaller than their values 
in the bulk solutions, that is, 4.7 and ~ 5 [28]. This significant decrease 
of the peak values for the Cs-0 aind Cs-H radial distribution functions is 
a consequence of the contact adsorption of the cation onto the Cm ion as 
discussed above in connection with the RDFs plotted in Figure 1. 

In Figure 2 the first peak of ^ciH is shifted by ~ 0.1 nm, about the 0-H 
bond length, compared with This suggests that the water molecules 

nearest to the anion have a preferred orientation with one 0-H bond pointing 
directly towards the anion. For the sodium cation, the first peak of guM is 
about 0.06 nm farther than the first peak of gnao- This shift corresponds to 
the configuration where one lone pair of the ST2 water points directly at the 
cation. 

These conclusions concerning the orientation of water molecules in the 
hydration shells of the ions are confirmed by the curves in Figure 3. Here we 
plot the average orientations of the water molecules as functions of the dis- 
twce between the oxygen atom of the water and the given ion. To mefisure 
the molecular orientation, we introduce the variable 6, which is the angle 
between the water dipole and the ion-0 connecting line. If one 0-H bond 
vector of the water molecule lies exactly on the line connecting the oxygen 
atom to the ion {0 — 125.25®) or if one lone pair vector points directly to¬ 
wards the ion = 54.75®), then the average < cos0 > would be ±0.58. 
As indicated by the first two curves (from the top) in Figure 3, the water 
molecules surrounding Cm fo be oriented with their dipoles pointing 
towards the center of the negatively charged sphere. Within the first hydra¬ 
tion shell for Na'*' and Cl”, the average < cos^ > remains approximately 
constant, ~ 0.65 for sodium cation and ~ —0.6 for the chloride anion. For 
purposes of comparison, we note that in a bulk solution simulation without 
the Cm molecular-ion, Heihzinger and Vogel [28] found that the mean values 
of co&$ are 0.69 for Na"*" and —0.65 for Cl". These values au-e slightly larger 
than the values obtained in this work. The orientational ordering of the water 
molecules surrounding the cesium ion is severely disturbed by the oppositely 
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charged C«o molecular ion. If the C«o ion were absent, these dipoles would 
have their orientations similar to th<»e of the water molecules around the 
Na'*’ cation (the third curve in Fig.3). 

4 Conclusions 

The first layer of water is oriented with its dipole pointing towards the center 
of the ’spherical’ doubly charged C«o anion. The cesium cation is in physical 
contact on the spherical surface about 10% of the time. The smaller ions 
occupy extensive diffuse layers around the Ceo molecular ion. No adsorption 
of the sodium ion onto the spherical surface is observed. Its distance of 
closest approach is one water molecule, with a separation of two or more 
water molecules more likely. The chloride ion appears to spend part of the 
time in a solvent-separated adsorption configuration and part of the time 
in the diffusion region. The hydration structures of the sodium cation and 
the chloride anion are not significantly perturbed by the presence of the 
Ceo molecule. The hydration shell of the cesium ion is incomplete and the 
orientations of the water molecules around Cs"*" are severely disturbed by the 
surface of Ceo- 
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Figure Captions 


Fig. 1 Radial distribution functions for water-C60 pair in the NaCl solution, 
water-C60 pair in the Cs solution, sodium-C60 pur, chloride-C60 pair, 
and cesium-CCO pair (from top to bottom). 

Fig. 2 Hydration structure of the ions. The dotted curves are for ion-0 
radial distribution functions and the solid curves are for ion-H radial 
distribution functions. 

Fig. 3 Average cos^ as a function of the distance between an ion and 
the oxygen atom. Curves from top to bottom represent the molecular 
orientations around C60 in the NaCl solution, around C60 in the Cs 
solution, Na+, Cl“, and Cs"*". 
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